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ABSTRACT 

Present-day massive galaxies are composed mostly of early-type objects. It is unknown 
whether this was also the case at higher redshifts. In a hierarchical assembling scenario 
the morphological content of the massive population is expected to change with time 
from disk-like objects in the early Universe to spheroid-like galaxies at present. In 
this paper we have probed this theoretical expectation by compiling a large sample of 
massive {Msteiiar 1^ lO^'^'^^To^Mo) galaxies in the rcdshift interval 0<z<3. Our sample 
of 1082 objects comprises 207 local galaxies selected from SDSS plus 875 objects 
observed with the HST belonging to the P0WIR/DEEP2 and GNS surveys. 639 of our 
objects have spectroscopic redshifts. Our morphological classification is done in the V- 
band restframe both quantitatively (using the Sersic index as a morphological proxy) 
and qualitative (by vistial inspection). Using both tcchniqties we find an enormous 
change on the dominant morphological class with cosmic time. The fraction of early- 
type galaxies among the massive galaxy population has changed from '~20-30% at z~3 
to ~70% at z=0. Elliptical galaxies have been the predominant morphological class 
for massive galaxies since only 2: ~ 1. 
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1 INTRODUCTION 

The present-day massive galaxy population is dominated 
by objects with early-type morphologies (e.g. Baldry et al. 
2004, Conselice et al. 2006). However, it is still unknown 
whether this was also the case at earlier cosmic epochs. Ad- 
dressing this question is key to our understanding of the 
physical processes that drive galaxy evolution, as galaxy 
morphology is directly linked to the evolutionary paths fol- 
lowed by these objects. In fact, a profound morphologi- 
cal transformation of the massive galaxy population is ex- 
pected within the currently most favoured galaxy formation 
scenario, the hierarchical model. For massive galaxies this 
model predicts a rapid formation phase at 2<z<6 domi- 
nated by a dissipational in-situ star formation fed by cold 
flows (Oser et al. 2010; Dekel ct al. 2009; Keres et al. 2005) 
and/or gas rich mergers (Ricciardelli et al. 2010; Wuyts et al. 
2010; Bournaud et al. 2011). At the end of this phase, mas- 
sive galaxies are expected to be more flattened and disk-like 

* E-mail: fernEindo.buitrago@roe.ac.uk 
t Scottish Universities Physics Alliance 



than their lower redsliift massive counterparts (Naab et al. 
2009). After this monolithic-like formation phase, massive 
galaxies are predicted to suffer a period of intense bombard- 
ment by minor satellites (Khochfar & Silk 2006; Hopkins et 
al. 2009; Feldmann et al. 2010; Oser et al. 2010) that may 
eventually transform the original disk-like population into 
the predominant present-day spheroid-like population. 

Although the above scenario is very suggestive of a deep 
morphological transformation of the massive galaxy popula- 
tion, there is no compelling observational evidence support- 
ing this scenario. However, some recent works suggests that 
this could be the case (e.g. van der Wei et al. 2011, Cameron 
et al. 2011, Wcinzirl ot al. 2011). Probing this transforma- 
tion is difficult from the observational point of view due to 
the scarce rmmbcr of massive galaxies at high-z. However, 
the advent of wide area and deep near infrared surveys (e.g. 
Dickinson et al. 2003, Scoville et al. 2007, Conselice et al. 
2011) have opened up the possibility of exploring a large 
number of these gala^xies up to high redshifts. In this paper 
we address, to the best of our knowledge for the first time, 
the issue of the morphological transformation of massive 
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galaxies using a statistical representative sample of nearly 
~1000 galaxies with Af* ^ 10^^ h^g Mq obtained from the 
SDSS DR7 (z~0; Abazajian et al. 2009), POWIR/DEEP2 
(0.2<z<2; Bundy et al. 2006, Conselice et al. 2007) and 
GNS (1.7<z<3; Conselice et al. 2011) surveys. We have al- 
ready conducted a morphological quantitative analysis of 
the above galaxies in previous papers (Trujillo et al. 2007; 
Buitrago et al. 2008) where we have provided clear evidence 
for a significant size evolution for these objects since z~3. 
However, a visual classification of these galaxies and an anal- 
ysis of their overall profile shape have been missing. In this 
paper we take advantage of the combined power of the vi- 
sual and quantitative morphological analysis to explore how 
the morphological content of the massive galaxy population 
has changed with redshift. 

The structure of the paper is as follows: Section [2] is 
devoted to the data description and its analysis. Section [3] 
presents our main results and in Section |4] we discuss them. 
At the end of this work we include an Appendix contain- 
ing the simulations we have performed to test the accuracy 
of our structural parameter determination in the GNS. On 
what follows, we adopt a cosmology with Qrn.=0.3, f2A=0.7 
and Ho=70 kms"^ Mpc~\ 



2 DATA 

To accomplish our objectives and to be statistically mean- 
ingful we need a large number of massive galaxies at all red- 
shifts within our study. Ideally we would also like to study 
all of our galaxies at a similar restframe wavelength range. 
This is the reason behind our choice of working with several 
different surveys. The imaging for the local Universe galaxy 
reference sample was obtained using the SDSS DR7 (Abaza- 
jian et al. 2009) although our sample was selected from the 
NYU Value-Added Galaxy Catalog (DR6). This catalog in- 
cludes sing le Sersic (1968) profile fits for 2.65 x lO" galaxies 
(Blanton et al. 2005), from which 1.1 x 10^ galaxies have 
spectroscopic information. Stellar masses come from Blan- 
ton & Roweis (2007), which uses a Chabrier (2003) IMF. 
We limited our work to all the massive (Af* ^ 10^^ h^^ Mq) 
galaxies with spectroscopic redshifts up to z = 0.03. We 
have selected this redshift as an upper limit as it contains 
a local sample with a number of objects (~200) similar to 
the number of galaxies we have in our higher redshift bins. 
On doing this we assure they are all affected statistically in 
a similar way. By selecting z=0.03 we also guarantee that 
these galaxies are retrieved from a sample that is complete 
in stellar mass. One object in this local sample was rejected 
£18 we discovered it was a stellar spike. Our final local sample 
contains 207 galaxies. We have used the g-band imaging of 
SDSS to classify visually our local sample. 

In the redshift range 0.2 < z < 2 we utilised the Palo- 
mar Observatory Wide-field InfraRed (POWIR)/DEEP2 
survey (Bundy et al. 2006, Conselice et al. 2007, 2008). For 
this part of the analysis, we restricted ourselves to the HST 
ACS I-band coverage in the Extended Groth Strip (ECS). 
The sample of massive galaxies selected from this survey 
constitutes the largest sample of massive galaxies in this 
redshift range published to date. The EGS field (63 Hubble 
Space Telescope tiles) was imaged with the Advanced Cam- 
era for Surveys (ACS) in the V(F606W, 2660s) and I-band 



(F814W, 2100s). Each tile was observed in 4 exposures that 
were combined to produce a pixel scale of 0.03 arcsec with a 
Point Spread Function (PSF) of 0.12 arcsec Full Width Half 
Maximum (FWHM). The depth reached is Iab = 27.52 (5cr) 
for point sources, and about 2 magnitudes brighter for ex- 
tended objects. Complementary photometry in the B, R and 
I bands was taken with the CFH12K camera at the CFHT 
3.6-m telescope and in the Ks and J bands with the WIRC 
camera at the Palomar 5-m telescope. Standard photomet- 
ric techniques (see Conselice et al. 2007) were performed in 
order to get both photometric redshifts and masses (using 
again Chabrier IMF). In total, 421 massive galaxies posses 
spectroscopic redshifts out of the total 795 in this survey. 
There were another 35 massive galaxies in the parent sam- 
ple, but they were excluded as they are identified as AGN 
and hence they may skew our results. 

For the highest redshift bins we used the GOODS NIC- 
MOS SurvejQ(GNS; Conselice et al. 2011). The GNS is a 
large HST NICMOS-3 camera program of 60 pointings cen- 
tered around massive galaxies at z — 1.7 — 3 at 3 orbits 
depth, for a total of 180 orbits in the F160W (H) band. 
Each tile (52"x52", 0.203" /pix) was observed in six expo- 
sures that were combined to produce images with a pixel 
scale of 0.1 arcsec, and a PSF of ~ 0.3 arcsec FWHM. Pho- 
tometric redshift and stellar masses were derived in a similar 
way as in POWIR. Some (11) spectroscopic redshifts were 
added a posteriori to the data (Popesso et al. 2008, Barger 
et al. 2009). This sample is the largest massive galaxies com- 
pendium (80 objects) at 1.7 < 2 < 3 we are aware of. 

2.1 Quantitative and visual morphological 
classification 

Once we selected the final sample of objects, the surface 
brightness distributions of all our galaxies were examined 
using GALFIT (Peng et al. 2010). We fit aU our objects with 
a single Sersic model convolved with the PSF of the images, 
selected using real stars. Difl^erent, non saturated, stars were 
used to check the reliability of our quantitative morphologi- 
cal analysis results. The exact details of our analysis can be 
found in Trujillo et al. (2007; POWIR sample) and Buitrago 
et al. (2008; GNS sample). In relation to the SDSS imaging, 
although the NYU catalog already provides structural pa- 
rameters obtained using Sersic fits to the galaxies, for the 
sake of consistency with our methodology, we ran GALFIT 
on the SDSS images of these galaxies to obtain structural 
parameters. It is known that the NYU catalog has a system- 
atical underestimation of the Sersic index, effective radius 
and total flux, as it is reported in the simulations performed 
in Blanton et al. (2005) and in the appendix of Guo et al. 
(2009). Our findings agree with this fact, as we find an offset 
of 26 ± 2% for the circularised effective radius values of our 
galaxies and another 14 ± 3% for the Sersic indices. 

In addition to the quantitative morphological analysis 
explained above, visual morphological classifications were 
derived for all the galaxies in our sample. To assure a high 
reliability in our results, two authors of this paper (FB and 
IT, with checks by CC) classified visually all the galaxies 
in an independent way. We divided our sample according to 
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Figure 1. Some examples illustrating our morphological criteria (columns) for different galaxies in our sample. Each row shows galaxies 
of the different surveys. Please note the different scales of each image due to each galaxy redshift (lower left corner); according to the 
cosmology used in this work, 10 arcsec in SDSS are ~ 6 kpc at z~0.03, while 1 arcsec in the HST imaging at 2; ^ 1 is ~ 8 kpc. Despite the 
decrease in angular resolution and the cosmological surface brightness dimming with redshift, the exquisite HST depth and resolution 
(~10 times better than ground-based SDSS imaging) allow us to explore the morphological nature of the high-z galaxies. Note that 
irregulars and mergers are in the same morphological class (peculiars). 



the Hubble classification scheme into spheroid-like objects 
(E-I-SO or early-type), disk-like objects (S or late-type) and 
peculiar galaxies (either irregular galaxies or ongoing merg- 
ers). In Figure[l]we show some examples of our classification 
scheme at different redshifts. Very conspicuous bulge sys- 
tems were identified as early-type objects. Both E and SO 
galaxies are hence included together in the same morpho- 
logical class. We avoid segregating between E and SO since, 
at high-z, it is a difficult task to distinguish between these 
two types of galaxies, and we prefer to remove this potential 
source of error. Spiral or late-type morphologies are detected 
by a central brightness condensation located at the centre 
of a thin disk containing more or less visible spiral arms 
of enhanced luminosity. Lastly we joined irregular (unsym- 
metrical) galaxies and mergers in the same class, again to 
avoid any missclassification at high-z where the differences 
between these types are more difficult to interpret. 

It is not straightforward to asses the robustness of our 
visual classifications as this is a pattern recognition problem 
and as such cannot be addressed by standard algorithm pro- 
cedures. Nevertheless, at z~0, we can compare our results 



with other alternatives coming from independent studies. 
First, we compare our results with the SDSS Bayesian auto- 
mated morphological classification by Huertas-Company et 
al. (2011). There are 190 out of our 207 galaxies in common 
where can make a direct comparison. They applied support 
vector machine techniques (Huertas-Company et al. 2008) 
to associate a probability to each galaxy being E, SO, Sab 
or Scd. For those galaxies where they have assigned a prob- 
ability larger than 90% of pertaining to a given class, their 
neural network agrees with our visual classification for 89% 
of the early-types and 68% of the late-types. Moreover, all 
our SDSS local galaxies have been visually classified within 
the Galaxy Zoo project (Lintott et al. 2011). We find that 
112 out of the 121 galaxies that we classified as early-type 
are classified as ellipticals by Galaxy Zoo (i.e. ~93%). For 
spiral galaxies we get 48 out of 62 (i.e. ~77%). The discrep- 
ancies in our late-type galaxies identifications arise from the 
difficulty to interpret SO galaxies without the astronomer's 
"trained-eye" intervention and, as stated previously, they 
are included as early-type objects throughout our study. 
Consequently, our local classification seems to be robust. 
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2.2 Potential observational biases 

We acknowledge, however, that at higher redshifts visual 
morphological classification is more controversial for several 
reasons. First, the cosmological surface brightness dimming 
may afi^ect the recognition of fainter galactic features and 
second, the angular resolution is poorer at higher redshift. 
Nonetheless, the first effect is compensated by the increase of 
the intrinsic surface brightness (e.g. Marchesini et al. 2007; 
Prescott, Baldry & James 2009, Cirasuolo et al. 2010) of the 
galaxies due to having a higher star formation rate in the 
past, and the fact that their stellar populations are younger. 
In relation to the angular resolution, at z=0.03, one arcsec 
is equivalent to ~0.6 kpc, whereas at l<z<3 it is ~8.0 kpc. 
Fortunatelly, the higher resolution imaging used for explor- 
ing the morphologies of our high-z galaxies (FH[WM~0. 1-0.3 
arcsec) compared to the local ones (FWHM~1.0-1.5 arcsec) 
alleviates this problem, although in general, a smoother sur- 
face brightness distribution due to the worse resolution is 
expected. All these effects combined imply that at higher 
redshifts there would be a larger number of featureless ob- 
jects that visually would be confused with early-type galax- 
ies. We will show in the next section that this is the opposite 
of what we find, ultimately giving stronger support to the 
results of this paper. 

2.3 Morphological K-corrections 

The K-correction effect is another potential source of er- 
ror both within the quantitative and visual morphological 
classification. We have selected our filters at each survey to 
minimise this effect, and to observe the galaxies as close as 
possible to the g-band restframe. Nonetheless, our classifica- 
tion at 1.3 ^ z ^ 2 could be compromised by using F814W 
as this filter is tracing the UV restframe of our targets. 

We explored how relevant this effect is by analysing the 
properties of 24 galaxies with z<2 in our POWIR/DEEP2 
sample which also have H-band NICMOS imaging. In Tru- 
jillo et al. (2007) we discussed the size difference between 
the optical and near-infrared for these galaxies (their Fig. 
4). We found a non-systematic bias, but a scatter of 32% for 
these measurements of the effective radius. With regards to 
the Sersic index, we found an offset of 30±9% towards larger 
indices in the H-band. The difference in the visual morphol- 
ogy between the I and the H-bands shows that 19 galaxies 
(79 ±8%) have the same morphology in the two filters, while 
only 5 (21 ± 8%) are catalogued differently. Our errors for 
visual classifications are the binomial standard deviations. 
We found in ACS 13 (54 ± 10%) early-type galaxies (10 in 
NICMOS, 42 ± 10%), 6 (25 ± 9%) late-type galaxies (10 in 
NICMOS, 42 ±10%) and 5 (21 ±8%) peculiar galaxies (4 in 
NICMOS, 17 ±8%). 

In addition to this analysis of galaxies in the EGS, we 
compared the difference between the I and H band mor- 
phologies for those galaxies in the GNS with 1.7<z<2 (which 
is the redshift interval where our POWIR/DEEP2 and GNS 
sample overlap). We used the I-band ACS imaging of the 
GOODS fields (Giavalisco et al. 2004). Postage stamp im- 
ages for 20 common galaxies were retrieved from the RAIN- 
BOW database (Barro et al. 2011). Our GALFIT analysis 
showed that the effective radius and the Sersic index are 
recovered without any significant offset, but with a large 



scatter as in the aforementioned Trujillo et al. (2007). Re- 
garding their visual morphologies, we found that 6 galax- 
ies (30 ± 10%) were not possible to classify reliably due 
to the few pixels that compose their image, most proba- 
bly due to dust obscuration (Buitrago et al. 2008, Bauer et 
al. 2011). For the remaining 14 galaxies, 11 (55 ± 11%) have 
the same visual morphology, while for 3 galaxies (15 ± 8%) 
there is a difference. For the detections in the ACS camera, 
4 (28 ± 12%) are early-types (6 for their NICMOS coun- 
terparts, 43 ± 11), 6 (43 ± 11%) are late-types (5 in NIC- 
MOS, 36 ± 13%) and 4 (28 ± 12%) are peculiars (3 in NIC- 
MOS, 21 ± 11%). Summarising, K-correction undoubtedly 
play a role, but visual morphologies are robust against these 
changes within our study. 

2.4 Axis ratios 

Finally, a further test we can conduct to quantify the ro- 
bustness of our visual classification is to explore the axis 
ratio distribution of our objects. The axis ratio distribu- 
tion of local disk galaxies has a mean value of ~0.5 (Ryden 
2004). On the other hand, the axis ratio distribution of the 
nearby E/SO population is known to peak at around 0.7-0.8 
(Ryden, Forbes & Terlevich et al. 2001). Figure [5] displays 
the distributions of axis ratios in our sample according to 
our visual morphological determinations. In Table [T] we also 
show the mean axis ratio for our different galaxy popula- 
tion as a function of redshift. We find that the objects that 
are visually classified as early-type galaxies have a typical 
axis ratio of ~0.7 (independent of their redshift). Also, for 
galaxies visually classified as disks, the axis ratio is indepen- 
dent of the redshift with an average b/a~0.55. Both values 
are in good agreement with the expectation from the local 
Universe. This test reinforces the idea that our visual clas- 
sifications are accurate. 



3 RESULTS 

The evolution of galaxy morphology with redshift can be ad- 
dressed in two different ways: quantitative (exploring how 
their structural parameters have changed with time), and 
qualitative (probing how the visual appearance has evolved 
with redshift). In the local Universe, the structural prop- 
erties of massive galaxies (mainly its light concentration) 
can be linked with their appearance. In particular, as a first 
approximation one can identify disk or late-type galaxies 
with those galaxies having lower values of the Sersic in- 
dex (n~l; Freeman 1970) and early-type galaxies with those 
having a profile resembling a de Vaucouleurs (1949) shape 
(n~4). This crude segregation based on the Sersic index 
was shown to work reasonably well by Ravindranath et al. 
(2004). Whether this equivalence also holds at higher red- 
shift is not clear, and in this paper we explore this issue. 

Table [l] lists the mean values of the structural parame- 
ters from the galaxies in our sample, splitting them among 
a number of redshift ranges and their visual morphologies. 
Concerning effective radius measurements, we retrieve again 
the reported size decrease for massive galaxies with redshift. 
As explained in Section [2.4[ axis ratios values per morpho- 
logical class are independent of redshift, which agree with 
our morphological classifications. Finally, there is a tendency 
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Figure 2. Axis ratio distributions for our sample of massive galaxies according to their visual morphology. Vertical dashed lines are the 
median values of each histogram (0.70, 0.54 and 0.59, respectively). While for disks the distribution is rather symmetrical, this is not 
the case for early-type objects, which peak at higher axis ratio values. Van der Wei et al. (2011) argued that massive galaxies with axis 
ratios ^ 0.5 are most probably related with late-type objects, as this histograms confirm. Interacting/peculiar galaxies should be taken 
out of these considerations. The histograms are also in agreement with the estimations of axis ratio distributions for local early- type 
and late-type objects (Ryden, Forbes & Terlevich et al. 2001; Ryden 2004). This fact highlight the reliability of our visual morphological 
classification. 



for the Sersic indices (observed using both mean and median 
values, see Fig.|3| to become progressively smaller at increas- 
ingly higher redshifts. Moreover, there is a separation in the 
average Sersic index for late and early-type massive galaxies 
at all redshifts. 

In Fig. |4] we show the Sersic index distribution for our 
different visually classified morphological types as a func- 
tion of redshift. At all redshifts, massive galaxies identified 
visually as late-types show low Sersic index values. This re- 
inforces the idea that the stellar mass density distributions 
of rotationally supported systems are close to an exponen- 
tial profile. However, the distribution of the Sersic index for 
these late-type galaxies show a tail towards larger values. 
This is normally interpreted as the result of the bulge com- 
ponent. In fact, the excess of fight caused by the bulge at 
the centre of the disk will increase the value of this concen- 
tration parameter when the galaxies are fitted just using a 
single Sersic model. Interestingly, we observe that at higher 
redshift the prominence of this tail of higher Sersic indices 
decreases for the late- type galaxies. 

One could be tempted to interpret this result as a con- 
sequence of the disappearance of prominent bulges at higher 
redshifts. However, a detailed exploration of this issue is be- 
yond the scope of this paper. In the same figure, we show 
the distribution of the Sersic index for massive galaxies vi- 
sually classified as early-types. We see that at low redshift, 
the distribution of Sersic indices for these galaxies predom- 
inantly show large values of concentration as it is expected. 
Up to z~1.5 there is a peak around n~4-6 (see also Table 
[T|. A general trend is also observed: there is a progressive 



shift towards lower and lower Sersic index values as redshift 
increases. 

We conducted a series of Kolmogorov-Smirnov (KS) 
tests in order to check the significance of the changes in 
the Sersic index with redshift. In Figure [7] we show the evo- 
lution of the KS significance levels for the different mor- 
phological types in our sample of massive galaxies. The KS 
test is a non-parametric method of comparing probability 
distributions. We used as a base comparison the lowest red- 
shift bin Sersic index distribution (left panel), and at the 
highest redshift one (right panel). Despite the uncertain- 
ties, mainly due to the number statistics, it seems that the 
distributions of Sersic indices gradually change from low to 
high redshift, and vice-versa. When we conduct the analy- 
sis using the high-z bin as the reference instead of the local 
one, the evolution of the Sersic index distribution is sta- 
tistically more significant on average. The fact that several 
redshift bins repeat significance level values in the left panel 
shows us the Sersic index distributions of early types are 
similar to the one at low-z, while this is not true for late- 
types. Possible explanations are the homology of early-type 
galaxies, and the double-peak which only exists in our local 
late-types Sersic indices, that is linked with the growth in 
the light profile wings. Consequently, the checks we carried 
out attempting to characterise the change of Sersic indices 
with redshift are positive (namely the Table [l] the Figure [s] 
and these KS tests), even though the KS comparison with 
the local Sersic index distribution is not as smooth as the 
others. The reason for this change or shift we observe could 
be either a real effect, produced by a decrease in the tail of 
the surface brightness distribution of the massive galaxies 
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Figure 3. Evolution of the mean Sersic index values over redshift according the visual classifications of the galaxies within our sample 
(see also Table Data points for early-types and peculiars are slightly offset for the sake of clarity. Dashed lines are similar but taking 
median values instead. Errors bars are the uncertainty of the mean (a/ ^ (N — 1), being a the standard deviation and N the total number 
of galaxies for each point). They are slightly larger at 0.2<z<0.6 because of the comparatively poor statistics at this redshift interval. 
Prom that epoch to higher redshifts there is a clear separation between early-type massive galaxies and the rest of visual types, being 
all the average Sersic indices lower at increasing redshift. 



at higher redshift, or an artificial one, produced by a bias at 
recovering large Sersic index values. 

To explore this last possibility we have conducted ex- 
tensive simulations to check whether there is any bias in 
the recovery of the Sersic index. In the case of the POWIR 
sample the simulations are fully explained in Trujillo et al. 
(2007). In this study we did not find any significant trend in 
either the sizes or the concentration of the galaxies (see their 
Fig. 3), except for a slight understimation of ~ 20% in the 
Sersic index of the very faint Iab > 24 spheroid-like galax- 
ies. We carried out a similar analysis here for the galaxies 
in the GNS sample. 

The results of these simulations are comprehensively ex- 
plained in the Appendix A at the end of this work. In sum- 
mary we find that for objects with disk-like surface bright- 
ness profiles (i.e ni„put<2.5), both sizes and Sersic indices 
are recovered with basically no bias down to our limiting 
explored H-band magnitude. However, by increasing the in- 
put Sersic index we find biases in the determination of the 
sizes and n. For a galaxy with ni„put~4 and H=22.5 mag 
(our typical magnitude within the GNS catalogue), the out- 
put effective radii are ~ 10% smaller and output Sersic in- 
dices are ~ 20% smaller than our input galaxies. The results 
of these simulations however show that the decrease in the 
Sersic index we observe from z~2.5 to z=0 for the spheroid- 
like population (which is around a factor of ~2) can not be 
explained fully as a result of the bias on recovering the Sersic 
index. 



We used the output of our simulations (Routput, le, output 
and rioutput) to estimate the intrinsic (input) values of our 
sample and provide a more accurate representation of the 
evolution of the Sersic indices at high-z. On average, the 
observed Sersic indices grow by ~10% after these changes, 
although the individual values per galaxy depend on its ex- 
act position in the 3D space defined by the output magni- 
tude, Sersic index and effective radius. In Figure [4] we im- 
plemented these corrections for the GNS, and also for the 
POWIR/DEEP2 data using Trujillo et al. (2007) simula- 
tions. For the sake of clarity. Figure [5] displays the same 
results but without any corrections in two highest redshift 
bins. Even after these corrections are applied the trend we 
observe towards lower Sersic indices at higher redshifts is 
maintained. In fact, the corrections are minor. We discuss 
and interpret the histograms of Fig.|4]in the next section. In 
relation to the distribution of the Sersic index for the galax- 
ies we classified as interacting or irregulars, we see a larger 
spread. 

Many studies (e.g. Shen et al. 2003; Barden et al. 2005; 
Mcintosh et al. 2005; Trujillo et al. 2006) have used n=2.5 as 
a quantitative way to segregate between early and late-type 
galaxies. We explore, using this criteria, how the percent- 
ages of the different types of massive galaxies evolve with 
redshift. This is shown in Fig. [6^. This figure clearly indi- 
cates that the fraction of massive galaxies with lower Sersic 
index values has dramatically increased at higher redshift. 
If the association between the Sersic index and the global 
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Figure 4. Sersic index distribution of massive (M* ^ 



) galaxies at different redshift intervals. Tlie Sersic indices of the 



individual galaxies have been corrected following the simulations presented in Trujillo et al. (2007; POWIR) and the Appendix of this 
paper (GNS). Color coding is related with visual morphology: blue for late-type galaxies, red for early-type galaxies and green for peculiar 
(irregulars/mergers) galaxies. It is pertinent to note that the sharp peak for late-type objects at 0.2<z<0.6 is due to the small number 
statitics for this class of massive galaxies at this redshift interval (see Table [ij. For our SDSS sample, the Sersic index of disky objects 
are mainly located between l<n<3 but for some galaxies extend up to n = 5. Conversely, the Sersic index of spheroid galaxies starts 
at n~3 and then peaks at n~5. The distributions of the Sersic index steadily shift to lower values at increasing redshift with the high 
Sersic index values progressively disappearing. 
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Figure 5. These are the highest redshift histograms of the Figure|4] showing the observed Sersic indices values, without any a posteriori 
correction based on Trujillo et al. (2007) or our current GNS simulations (Appendix A). The more noticeable change is seen for the GNS 
data, where it is very conspicuous the non-existence of any large (n > 4.62) Sersic index. The difference between these histograms and 
the ones presented in Figure]?] is small. 
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morphological type that holds at low redshift also applies at 
high-z this would imply that massive galaxies at the high- 
z universe were mostly late-type (disk) galaxies. However, 
there is no guarantee that such an association holds at all 
redshifts. For this reason, we explore the evolution of the 
fraction of different galaxy types with redshift using the vi- 
sual morphologies (see Fig. [gJj) . We find that the population 
of visually classified massive disk galaxies remains almost 
constant with redshift with perhaps a slight (if any) increase. 
The most dramatic changes are associated with the early- 
type and irregular/mergers classes. The fraction of visually 
classified E/SO galaxies increased by a factor of three since 
z~3 to now, whereas a reverse situation is seen for the irreg- 
ular/merging galaxies. This latter fact agrees with merging 
becoming more important in massive galaxy evolution at 
increasing redshift (Conselice et al. 2009, Bluck et al. 2009, 
2011). At z~2.5 late-type and peculiar objects account for 
the majority of massive galaxies. One of the most important 
outcome of Fig. [6] is that the E/SO type has been the dom- 
inant morphological fraction of massive galaxies only since 
z~l. 

The number density of massive galaxies has significantly 
changed since z~3 (e.g. Rudnick et al. 2003; Perez-Gonzalez 
et al. 2008, Mortlock et al. 2011, Conselice et al. 2011) with 
a continuous increase in the number of these objects in the 
last ~11 Gyr. In order to probe the emergence of the differ- 
ent galaxy types explored in this paper we have estimated 
the number density evolution of each class. To do this, we 
have used the Schechter fits to the stellar mass functions 
provided by Perez-Gonzalez et al. (2008; their Table 2). 
We have integrated these functions for all massive objects 
with Msteiiar 5s W^^h^QMQ. We have later multiplied those 
numbers by the fractions we have estimated for the differ- 
ent classes of galaxies explored in this work. We show these 
number density evolution in Fig.[6j; &[6jl. The number den- 
sity of both disk-like and spheroid-like massive galaxies, ac- 
cording to their Sersic index, have changed with time. This 
evolution is particularly significant for spheroid-like objects 
which are now a factor of ~10 more numerous per unit vol- 
ume than at z~2. The number of massive disks has also in- 
creased as cosmic time progresses, but at a lower rate than 
spheroid galaxies. Finally, the number density of massive ir- 
regular/merging galaxies has grown only very midly, if at 
all, in the last ~11 Gyr. 



4 DISCUSSION 

The evidence collected in the previous section suggests 
that there is a strong evolution in the morphological prop- 
erties (both quantitative and qualitative) of the massive 
galaxy population over time. At high redshift, in agreement 
with the theoretical expectation, the dominant morpholog- 
ical classes of galaxies are late-types and peculiars. Gonse- 
quently, the morphological and structural types of the ma- 
jority of the massive galaxies at a given epoch has dramat- 
ically evolved as cosmic time increases. Two effects could 
play a role towards explaining this significant change of the 
dominant morphological class. On one hand, the galaxies 
that are progressively been added into the family of mas- 
sive objects (i.e. by the merging of less massive galaxies) are 
incorporated with already spheroidal morphologies. On the 



other hand, the already old massive galaxies can also evolve 
towards spheroidal morphologies due to frequent mergers. 
For instance, frequent minor mergers (Lopez-Sanjuan et al. 
2010, Kaviraj 2010, Lopez-Sanjuan et al. 2011, Bluck et al. 
2011) with the massive galaxy population will destroy ex- 
isting stellar disks, and also would be responsible for the 
appearance of long tails in their luminosity profiles. This 
scenario could explain why the evolution towards spheroid- 
like morphologies is stronger when we use the Sersic index 
n instead of the visual classification. In fact, the surface 
brightness of nearby massive ellipticals are well described 
with large Sersic indices due to their bright tails or en- 
velopes. These wings, however, seem to disappear at higher 
and higher redshifts (see Table [l] or Figure [sj just leaving 
the inner (core) region of the massive galaxies (Bezanzon 
et al. 2009; Hopkins et al. 2009; van Dokkum et al. 2010; 
Garrasco, Gonselice & Trujillo 2010). The disappearance of 
these tails is also connected with the dramatic size evolu- 
tion reported in previous works (see e.g. Trujillo et al. 2007; 
Buitrago et al. 2008; Trujillo, Ferreras & De la Rosa 2011). 
Gonsequently, it is not only that the typical morphology of 
the massive galaxy population is changing with redshift, but 
also that there is a progressive build-up of their tails, making 
the morphological evolution appears more dramatic when we 
use the Sersic index instead of the visual classification as a 
morphological segregator. 

If we were just using the information contained in the 
change of the fraction of morphological types with redshift 
we would be tempted to explain the morphological evolution 
as being just a consequence of a transformation from one 
class to another, however, the evolution in the number den- 
sity of all the classes suggests a more complex scenario. In 
fact, one of the results we can conclude from the evolution of 
the number densities of all the massive galaxy classes is that 
high-z massive disk-like galaxies can not be the only pro- 
genitors of present-day massive spheroid-like galaxies. They 
are just simply not enough in number to explain the large 
increase of the number density of elliptical galaxies at low 
redshifts. 

All the morphological classes, with the maybe exception 
of irregular/merging galaxies, have increased their number 
densities with cosmic time. These irregular/merging galaxies 
also have shorter structural life-spans than the other types, 
and therefore the galaxies in this class must be replenished 
through time. This emergence of massive galaxies is more 
efficient (by a factor ~2) at creating spheroid-like galaxies 
than disk-like objects from z~l to now (see Fig. |6|. The 
reason why the formation of elliptical galaxies is more ef- 
ficient at recent times than it was in the past is possibly 
linked to the less availability of gas during mergers creat- 
ing new galaxies (Khochfar & Silk 2006, 2009, Shankar et 
al. 2010, Eliche-Moral et al. 2010), and thus more likely to 
contain denser and more concentrated light profiles (Bluck 
et al. 2011). 



5 SUMMARY 

Using a large compilation of massive (M lO^^/iyQ^M©) 
galaxies (~1100 objects) since z~3 we have addressed the 
issue of the morphological change of this population with 
time. We have found that there is a profound transformation 
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Figure 6. Panel A): Fraction of massive {Af» ^ lO^^h^^MQ) galaxies showing disk-like surface brightness profiles (n < 2.5) and 
spheroid-like ones (n > 2.5) as a function of redshift. Different color backgrounds indicate the redshift range expanded for each survey: 
SDSS, POWIR/DEEP2 and GNS. Error bars are estimated following a binomial distribution. Sersic indices are corrected based on our 
simulations (Trujillo et al. 2007 and Appendix A in the present paper). B): Same as Panel A) but segregating the massive galaxies 
according to their visual morphological classification. Blue color represents late type (S) objects and red early type (E-I-SO) galaxies, 
while peculiar (ongoing mergers and irregulars) galaxies are tagged in green. Panel C): Number density evolution of massive galaxies 
splitted depending on the Sersic index value. The solid black line corresponds to the sum of the different components. Panel D): Same 
as panel C) but segregating the massive galaxies according to their visual morphological type. 



in the morphological content of massive galaxies during this 
cosmic interval. Massive galaxies were typically disk-like in 
shape at z 1 and elliptical galaxies have been only the 
predominant massive class since that epoch. The fraction 
of early-type morphologies in massive galaxies has changed 
from ~20-30% at z~3 to ~70% at z=0. 

We have addressed the morphological transformation of 
the massive galaxies using a quantitative approach, based 
on GALFIT fits to the surface brightness distribution of the 
galaxies, and a qualitative approach based on visual classifi- 
cations. Both analysis agree on a clear morphological change 
in the dominant morphological class with time. In particu- 



lar, the quantitative approach, which uses the Sersic index 
as a morphological segregator, shows that the number of 
galaxies with low Sersic index at high-z was higher than in 
the present day universe. We interpret this as a consequence 
of two phenomena: a decrease in the number of early-type 
galaxies at higher redshift plus an intrinsic decrease of the 
Sersic index values of those elliptical galaxies at earlier cos- 
mic times due to the lost of their extended envelopes. 
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Redshift Redshift 

Figure 7. Kolmogorov-Smirnov significance levels comparing the SDSS local sample Sersic index distribution with all the other redshift 
bins (left panel) and the GNS highest redshift bin Sersic index distribution against the rest of the distributions (right panel). The color 
coding is the same as in the previous figures: red is for early-type galaxies, blue is for late-type galaxies and green is for peculiar galaxies. 
Dashed lines represent the same results but taking into account the corrections for Sersic indices inferred from our simulations. The 
reason why dashed lines encompass the whole plot in the right panel is that the last distribution (GNS) is available with and without 
corrections, but this does not occur for the local SDSS sample. The general trends show that, within our uncertainties, the distributions 
are diverging from the local relation (left panel) or progressively converging up to the highest redshift bin (right panel). The distributions 
with the Sersic indices corrected according to our simulations show a lower departure with respect to the fiducial distribution, but still 
a significant one. 
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Early-type galaxies 



Redshift Range 
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APPENDIX A: GOODS NICMOS SURVEY 
MASSIVE GALAXIES SIMULATIONS 

The purpose of this Appendix is to explore the robustness of 
the structural parameters of the massive galaxies (1.7<z<3) 
in our GNS sample (H-band, F160W filter, HST NICMOS-3 
camera, 3 orbits depth; Conselice et al. 2011). As explained 
in the main text of the paper, a set of simulations similar to 
the ones presented here were already conducted for the ACS 
imaging used to analyze the galaxies in the redshift inter- 
val 0.2<z<2 (Trujillo et al. 2007). To identity the ranges of 
the structural parameters to explore in our simulations, we 
use as a guide the ranges found in the quantitative morpho- 
logical analysis based on GALFIT of the real GNS massive 
galaxies (Buitrago et al. 2008). These were: 

0.15 < re{arcsec) < 0.61 
0.34 < n < 4.62 
0.19 <ar < 0.92 
-84.03 <pa< 85.28 

20.5 < HAB{mag) < 24 

where r^, n, ar, pa, Hab stand for effective radius, Sersic 
index, axis ratio, position angle and derived i/As-band mag- 
nitude. The only exceptions were one galaxy with n = 0.17 
and two others with 24 < Hab < 24.5. Taking these quan- 
tities into account, we simulated 16000 galaxies utilising the 
IDL routines built for Haufiler et al. (2007), with the struc- 
tural parameters randomized within these ranges: 

0.15 < re{arcsec) < 2.0 
0.25 < n< 8.0 
0.1 <ar < 1.0 
-89.99 <pa< 90.0 

20 < HAB{mag) < 25 

The structural parameters of the mock galaxies were dis- 
tributed linearly along the full parameter space, except for 
the effective radii which were logarithmically sampled as we 
specially wanted to explore objects with small angular radii 
due to the observed compactness of massive galaxies at high- 
z. 

Images of every mock galaxy were created placing these 
objects randomly in the GNS pointings. We placed a mock 
galaxy on each GNS poiting for every simulation in order to 
avoid altering the typical density (i.e. number of neighbour 
galaxies) of the GNS imaging. Each model galaxy (i.e. the 
2D surface brightness distribution following the Sersic func- 
tion) was convolved with a representative PSF. Specifically 
we used one of the five natural stars which were utilised in 
Buitrago et al. (2008). To obtain errors in the same way as 
in that paper, we also ran GALFIT using these five different 
stars and then taking the mean values. 

We have identified that the main source of uncertainty 
in the NICMOS data is the change of the PSF among the 
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different pointings. To illustrate how this affects the recovery 
of the structural parameters we first present in Figure Hjand 
[9] how our parameters are recovered when we use the same 
PSF for creating and recovering the mock galaxies. In Figure 
[T0land[TT| we show the effects on the parameters when we 
compare the input values with the average values obtained 
using the five different PSFs. 

In Figure [To] we show the relationship between the rela- 
tive errors in the structural parameters (magnitude, effective 
radius and Sersic index) versus the galaxy input magnitude. 
The relative errors are calculated as (output-input) /mput (i. 
e., a negative % refer to cases where the output is smaller 
than the input and vice- versa). The left column of the plot 
displays the structural parameters of individual mock galax- 
ies, whereas the right column shows their means in bins of 
0.5 mag. The mean values of the structural parameters were 
derived using a robust method which removes the 5cr out- 
liers. Error bars represent the standard deviation of the sam- 
ple. To appreciate how the effect on the structural parameter 
is linked to the input Sersic index of the mock galaxies we 
split the sample into four groups (0 < n < 2, 2 < n < 4, 



4 < n < 6 and 6 < n < 8). The results shown in Fig. 10 are 
tabulated in the Table |4] At increasing Sersic index values, 
the recovery of the structural parameters is more affected. 
We note that galaxies with low Sersic index are well recov- 
ered down to our faintest magnitude. An average galaxy in 
our GNS sample (H=22.5 mag) with an input Sersic index 
of n=4 will have its effective radius biased only by a ~15% 
and its Sersic index around ^ 25%. 

In addition to the dependence of the apparent magni- 
tude of the objects on recovering their structural parame- 
ters, in Figure [TT] we explore what is the effect of the size 
(lower row) and intrinsic shape (upper row) for this matter. 
Galaxies are colour coded in these figures according to their 
magnitude. Combining the information contained in Fig. |10| 
and |ll| we find that the key parameters for retrieving accu- 
rate structural parameters are the apparent magnitude and 
the Sersic index. The effective radius of the objects play a 
minor role. The results show on Fig. [Tl] are tabulated in the 
Tabled 

The most interested output of our tests using the mean 
value from the fits of different PSFs is that the size of the 
source plays now a fundamental role at characterising the er- 
ror on the structural parameters. As expected, large sources 
are less affected by changes in the PSF and the bias on the 
structural parameters remain basically the same as when 
we use just a single PSFs. However, at smaller sizes the 
effect of not knowing accurately the PSF at the source im- 
plies that the Sersic index uncertainty is large, although sizes 
are retrieved accurately. Summarising, neither any effect or 
combination of effects is large enough to modify the main 
results of this paper. Moreover, as stated on the main text 
of this paper, we use these simulations to correct, based on 
the observed (output) apparent magnitude, effective radius 
and Sersic index, the structural parameters presented in this 
work. 
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Table 2. Relative errors (%) and standard deviations on the structural parameters depending on the apparent magnitude (see Fig. 





< n < 2 


2 < n < 4 


4 < n < 6 


6 < n < 8 


20.0 < HAB,znput{mag) < 20.5 


SL/L 
Sre/r^ 
5n/n 


0.39 ± 1.75 
0.09 ± 2.03 
-0.34 ±4.08 


-0.53 ±3.21 
-1.16 ± 5.27 
-3.22 ±6.20 


-1.85 ±6.24 
-4.63 ± 12.35 
-6.25 ± 10.59 


-2.62 ± 7.36 
-7.19 ± 15.54 
-7.62 ± 12.07 


20.5 < HAB,r^put{mag) < 21.0 


5L/L 

Sre/re 
5n/n 


0.21 ±2.61 
0.01 ± 2.98 
-0.97 ± 5.63 


-0.99 ± 5.46 
-2.33 ±8.31 
-3.51 ± 8.48 


-2.95 ± 7.61 
-6.78 ± 14.35 
-7.69 ± 11.25 


-3.43 ±8.20 
-9.16 ± 17.30 
-9.13 ± 12.68 


21.0 < HAB,znput{mag) < 21.5 


SL/L 
Sre/r^ 
5n/n 


-0.12 ±3.85 
-0.40 ±4.58 
-1.42 ±8.21 


-1.34 ± 6.93 
-2.81 ± 11.39 
-5.65 ± 12.94 


-2.60 ± 10.14 
-5.08 ±20.67 
-7.72 ± 15.04 


-3.62 ± 10.74 
-7.85 ±21.32 
-9.08 ±16.12 


21.5 < H AB ,input(mag) < 22.0 


SL/L 

SVe/re 

5n/n 


-0.05 ±5.17 
-0.50 ±4.79 
-2.47 ± 9.68 


-2.56 ± 10.93 
-4.98 ± 17.21 
-6.85 ± 14.97 


-3.76 ± 13.39 
-8.46 ±25.80 
-10.83 ± 20.77 


-4.95 ± 13.55 
-9.54 ±29.76 
-11.80 ±21.14 


22.0 < HAB,znput{mag) < 22.5 


SL/L 

Sre/Ve 

5n/n 


-1.03 ± 7.28 
-1.70 ±8.48 
-3.55 ± 17.15 


-2.87 ± 12.06 
-5.33 ± 17.56 
-8.50 ± 19.65 


-5.83 ± 15.87 
-8.22 ± 27.02 
-11.86 ± 22.91 


-6.79 ± 18.50 
-14.28 ±32.00 
-17.65 ± 26.70 


22.5 < HAB,^nput(mag) < 23.0 


SL/L 
Sre/r-a 
5n/n 


-1.46 ±9.27 
-1.48 ± 10.75 
-3.74 ± 22.80 


-3.40 ± 18.16 
-5.55 ± 26.98 
-8.20 ± 25.31 


-6.97 ± 19.19 
-10.61 ± 33.02 
-13.99 ± 28.88 


-11.04 ±22.58 
-18.60 ±37.86 
-21.90 ±29.83 


23.0 < HAB,input{'mag) < 23.5 


SL/L 
Sre/r^ 
Sn/n 


-2.70 ± 18.92 
-3.38 ± 18.65 
-5.29 ± 29.66 


-5.54 ± 22.91 
-7.10 ± 33.57 
-11.14 ± 34.34 


-10.57 ± 22.24 
-15.11 ± 38.94 
-18.92 ±32.87 


-14.38 ±25.00 
-24.55 ± 39.17 
-29.79 ±31.89 


23.5 < HAB,input{mag) < 24.0 


SL/L 
Sre/r^ 
Sn/n 


-1.28 ±21.90 
-3.81 ± 26.00 
0.61 ± 37.18 


-8.12 ±22.41 
-10.21 ± 33.56 
-18.35 ±36.38 


-11.57 ± 26.88 
-17.86 ±43.47 
-24.90 ± 37.94 


-17.93 ±26.23 
-32.17 ±39.63 
-35.02 ± 33.59 


24.0 < HAB,m,put{mag) < 24.5 


SL/L 
Sre/re 
Sn/n 


-0.99 ± 28.30 
-2.13 ±36.39 
-7.57 ± 43.98 


-7.27 ± 34.80 
-13.48 ±41.48 
-27.35 ±41.82 


-16.06 ±33.78 
-29.16 ±44.22 
-39.53 ±40.77 


-15.98 ±34.87 
-32.70 ±42.43 
-44.10 ±36.75 


24.5 < HAB,r-nput{mag) < 25.0 


SL/L 

SVe/re 

Sn/n 


12.20 ± 51.04 
-12.85 ±45.18 
-11.73 ±49.37 


2.50 ±51.58 
-22.74 ± 46.45 
-38.99 ±42.35 


-5.73 ±44.09 
-31.77 ± 43.55 
-44.67 ± 42.69 


-15.50 ±43.68 
-36.15 ±48.64 
-47.58 ±40.68 
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Table 3. Relative errors (%) and standard deviations on the structural parameters (see Fig.^ 



Sr^/re 


20 < HAB(mag) < 21.5 


21.5 < HAsimag) < 23 


23 < HAB{mag) < 25 


< n < 2 
2 < n < 4 
4 < n < 6 
6 < n < 8 


-0.04 ±3.09 
-2.14 ± 8.27 
-5.71 ± 15.87 
-8.16 ± 18.05 


-1.00 ± 7.97 
-5.13 ±20.61 
-9.08 ±28.73 
-13.98 ±33.45 


-4.89 ±31.91 
-12.66 ±38.79 
-22.12 ±42.84 
-31.04 ± 42.55 


0.15 < re{arcsec) < 0.3 
0.3 < re{arcsec) < 0.6 
0.6 < re{arcsec) < 0.9 
0.9 < r^{arcsec) < 2.0 


-0.94 ± 5.45 
-2.28 ±9.02 
-4.15 ± 12.67 
-8.44 ± 18.68 


0.05 ± 13.79 
-2.71 ± 20.11 
-9.55 ± 24.26 
-15.66 ±32.69 


-0.87 ±32.27 
-10.47 ± 35.21 
-17.50 ± 40.12 
-35.02 ±43.53 


Sn/n 


20 < Hab < 21.5 


21.5 < Hab < 23 


23 < Hab < 25 


< n < 2 
2 < n < 4 
4 < n < 6 
6 < n < 8 


-0.84 ±6.23 
-3.75 ±8.71 
-7.08 ± 12.33 
-8.85 ± 14.25 


-3.57 ± 16.86 
-7.83 ± 20.63 
-12.19 ± 24.38 
-16.94 ± 26.32 


-5.51 ± 39.42 
-22.65 ±39.74 
-30.01 ± 39.47 
-38.55 ± 36.30 


0.15 < re{arcsec) < 0.3 
0.3 < re{arcsec) < 0.6 
0.6 < re{arcsec) < 0.9 
0.9 < reiarcsec) < 2.0 


-5.99 ± 12.04 
-3.93 ± 10.03 
-5.08 ± 11.05 
-6.38 ± 12.51 


-7.56 ±22.54 
-7.20 ±21.14 
-11.31 ± 21.56 
-13.74 ± 25.47 


-17.43 ± 39.18 
-19.31 ± 38.89 
-26.36 ±39.79 
-32.49 ±41.78 
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Table 4. Relative errors {%) and standard deviations on the structural parameters depending on the apparent magnitude using five 
different PSFs (see Fig.jTo}. 





< n < 2 


2 < n < 4 


4 < n < 6 


6 < n < 8 


20.0 < HAB,znput{mag) < 20.5 


5L/L 

Sre/re 
5n/n 


-0.15 ±2.95 
0.79 ± 5.25 
-4.39 ± 10.97 


-3.39 ±5.48 
-4.95 ±9.28 
-13.84 ± 16.60 


-9.09 ± 7.90 
-14.69 ± 14.79 
-24.76 ± 19.20 


-14.50 ± 8.47 
-26.56 ± 16.21 
-36.22 ± 19.58 


20.5 < HAB,input{mag) < 21.0 


5L/L 
Sre/r^ 
5n/n 


-0.14 ±3.82 
0.20 ± 5.36 
-4.16 ± 11.25 


-3.44 ± 6.32 
-4.57 ± 10.89 
-13.82 ± 17.91 


-9.31 ± 9.39 
-16.51 ± 16.28 
-23.84 ± 19.20 


-15.08 ±9.23 
-26.97 ± 17.87 
-38.19 ± 19.64 


21.0 < HAB,znput{mag) < 21.5 


SL/L 
Sre/r^ 
5n/n 


-0.57 ± 4.92 
-0.86 ± 7.00 
-5.76 ± 12.02 


-3.68 ± 7.55 
-5.73 ± 13.66 
-15.86 ± 17.54 


-9.49 ±9.25 
-16.66 ± 18.67 
-26.63 ± 20.50 


-14.29 ± 10.06 
-26.05 ± 18.35 
-36.08 ± 19.52 


21.5 < HAB,rnput(mag) < 22.0 


5L/L 
Sre/re 
5n/n 


-0.49 ±5.53 
-0.29 ± 7.33 
-6.04 ± 13.43 


-4.59 ± 10.36 
-7.05 ± 17.72 
-16.61 ± 18.69 


-9.94 ± 12.10 
-17.91 ± 22.02 
-27.84 ± 22.26 


-14.68 ± 11.82 
-25.85 ± 24.09 
-36.30 ±22.71 


22.0 < HAB,irtpudmag) < 22.5 


5L/L 
Sre/r^ 
5n/n 


-1.47 ± 7.60 
-1.75 ±9.51 
-8.18 ± 17.99 


-4.71 ± 13.41 
-7.38 ± 19.43 
-18.82 ± 21.73 


-11.55 ± 13.82 
-16.97 ± 25.13 
-28.87 ± 22.24 


-14.61 ± 15.32 
-26.68 ±26.94 
-39.18 ±24.90 


22.5 < H AB ,input(mag) < 23.0 


&L/L 
Sn/n 


-1.67 ± 10.03 
-1.28 ± 11.40 
-8.25 ± 22.92 


-5.89 ± 16.50 
-9.24 ± 26.02 
-18.29 ± 26.96 


-11.89 ± 18.38 
-17.99 ± 29.22 
-29.25 ± 27.00 


-17.96 ±21.49 
-27.32 ± 34.57 
-41.66 ±26.64 


23.0 < HAB,input{mag) < 23.5 


5L/L 
5n/n 


-3.80 ± 19.13 
-2.84 ± 21.33 
-7.49 ± 29.62 


-7.12 ± 22.24 
-10.78 ±32.29 
-20.09 ± 33.28 


-14.63 ± 19.75 
-20.79 ±34.39 
-32.47 ± 31.11 


-16.61 ± 24.12 
-25.78 ±39.19 
-42.58 ±29.90 


23.5 < HAB,^nput{mag) < 24.0 


SL/L 
Sre/r^ 
Sn/n 


-1.89 ± 22.85 
-2.21 ± 27.87 
-3.87 ± 38.56 


-10.19 ± 22.94 
-11.57 ± 35.12 
-24.95 ±36.87 


-14.19 ± 25.32 
-20.76 ±39.58 
-36.53 ± 34.02 


-20.52 ± 24.32 
-31.83 ±41.51 
-46.46 ±31.09 


24.0 < HAB,r-nput{mag) < 24.5 


SL/L 

SVe/re 

Sn/n 


-1.50 ± 25.40 
-3.60 ±35.90 
-10.32 ±42.23 


-10.23 ±32.27 
-14.78 ±43.22 
-32.03 ±42.13 


-16. 10± 34.15 
-30.37 ± 42.96 
-47.21 ± 38.40 


-16.71 ± 36.33 
-30.21 ±45.39 
-51.98 ±36.24 


24.5 < HAB,ir,put(mag) < 25.0 


SL/L 
Sre/r^ 
Sn/n 


3.48 ±38.67 
-11.83 ±45.10 
-17.46 ±48.50 


-0.99 ±47.79 
-24.02 ±44.12 
-42.41 ± 46.73 


-7.13 ±48.32 
-33.58 ±47.25 
-45.35 ± 44.40 


-17.63 ±39.12 
-36.64 ± 47.00 
-53.14 ±39.50 
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Table 5. Relative errors (%) and standard deviations on the structural parameters using five different PSFs (see Fig. [llj 



Sr^/re 


20 < HAB{mag) < 21.5 


21.5 < HAsimag) < 23 


23 < HAB{mag) < 25 


< n < 2 
2 < n < 4 
4 < n < 6 
6 < n < 8 


0.06 ± 5.77 
-5.12 ± 11.07 
-15.72 ± 16.22 
-26.47 ± 17.73 


-1.16 ±9.59 
-8.03 ±21.95 
-17.63 ± 25.57 
-26.40 ± 29.20 


-4.66 ±32.39 
-14.81 ± 38.86 
-25.64 ± 41.12 
-30.86 ± 43.40 


0.15 < re{arcsec) < 0.3 
0.3 < r^iarcsec) < 0.6 
0.6 < r^iarcsec) < 0.9 
0.9 < rs{arcsec) < 2.0 


-4.88 ± 13.91 
-12.01 ± 15.03 
-14.19 ± 17.17 
-17.15±21.19 


-2.74 ± 16.53 
-10.85 ± 19.95 
-16.98 ± 24.64 
-21.85 ±30.48 


-0.63 ±31.95 
-12.85 ±34.54 
-19.68 ±40.32 
-37.75 ±43.51 


Sn/n 


20 < Hab < 21.5 


21.5 < Has < 23 


23 < Hab < 25 


< n < 2 
2 < n < 4 
4 < n < 6 
6 < n < 8 


-4.89 ± 11.68 
-14.57 ± 17.28 
-25.06 ± 19.75 
-36.86 ± 19.70 


-7.60 ± 18.65 
-17.79 ± 22.82 
-28.64 ± 23.92 
-39.00 ± 24.87 


-9.07 ±39.19 
-29.09 ±40.43 
-39.53 ±37.22 
-48.32 ± 34.49 


0.15 < re{arcsec) < 0.3 
0.3 < re{arcsec) < 0.6 
0.6 < re{arcsec) < 0.9 
0.9 < reiarcsec) < 2.0 


-37.61 ± 22.17 
-21.64 ±20.35 
-16.58 ± 18.55 
-13.81 ± 17.93 


-37.41 ± 24.34 
-22.08 ± 24.19 
-19.83 ± 23.63 
-19.11 ± 25.80 


-35.55 ±39.30 
-29.12 ±39.11 
-29.86 ±40.01 
-35.62 ±42.06 
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Figure 8. Relative errors - (output-input) /input - of the structural parameters (magnitude, effective radius and Sersic index) of our 
simulated GNS galaxies. The right column shows the means in bins of 0.5 mag (with a 5(t outlier-resistant determination), being the 
error bars the standard deviation of the sample. The information in this plot is tabulated in Table |2] 
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Figure 9. Relative errors - output-input/input - of the effective radius (first column) and the Sersic index (second column) as a function 
of the input Sersic index (first row) and the input effective radius (second row). Galaxies are coloured according to their magnitude. 
For the sake of clarity, mean values (derived with a 5a outlier-resistant determination) where added using 4 intervals in effective radius 
and Sersic index, with the error bars being their standard deviation. Sersic index intervals are 0<n<2, 2<n<4, 4<ri<6 and 
6 < n< 8. Effective radius intervals are 0.15" < < 0.3", 0.3" < Ve < 0.6", 0.6" < re < 0.9" and 0.9" < < 2". Note that the colour 
of these mean points is the same as the one of the galaxy individual points. The information in this plot is tabulated in Table [s] 
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Figure 10. Same as in Figure|8] but using this time as the output parameters the mean values of the fits retrieved based on 5 different 
natural stars as PSFs. The results of this figure are tabulated in Table [4] 
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Figure 11. Same as Figure [9] but using this time as the output parameters the mean values of the fits retrieved based on 5 different 
natural stars as PSFs. The results of this figure are tabulated in Table [s] 



